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Abstract. We report on the BeppoSAX observation of 
the X-ray source X 1724—308 in the globular cluster 
Terzan 2. The broadband spectrum can be described as 
the superposition of a power-law with photon spectral in- 
dex r ~ 1.6 — 1.9 and a thermal component with typical 
temperature kT -1 keV. X 1724-308 is detected in the 
PDS up to 150 keV with a S/N ratio >3, notwithstand- 
ing a sharp exponential cutoff with E cu toff — 90 keV. The 
broadband spectrum can be best interpreted as due to 
the Comptonization of a Wien photon distribution with 
kTw ^T keV by a spherically-symmetric hot plasma with 
electron temperature kT e ~ 30 keV and optical depth 
r ~ 3. A comparison with the fainter and steeper state, 
which was observed by EXOSAT, suggests a correlation 
between the temperature of the Comptonizing electron 
distribution (or the plasma optical depth) and the X-ray 
intensity. The X-ray thermal component could originate 
from a boundary layer between the accretion disk and the 
compact object, or from the accretion disk itself. 
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1. Introduction 

The globular cluster Terzan 2 has been known to contain a 
luminous X-ray source since the late 1970s. An X-ray burst 
from a region including the cluster was observed by OSO-8 
(Swank et al. 1977), and the same mission detected persis- 
tent emission from a location consistent with the cluster. 
Uhuru (Forman et al. 1978) and Einstein High Resolution 
Imager (HRI, Grindlay et al. 1980) observations placed 
this source in the core of the cluster. The source is known 
either as X 1724-308 (Parmar et al. 1989; hereafter P89) 
or IE 1724—305. It has a rather hard (photon index T —2) 
power-law spectrum, with a 1-20 keV flux which varies by 
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about a factor of 4 (2.5 - 10 x 10~ 10 erg s _1 cm -2 ). Ein- 
stein data suggest that in a higher flux state (i 7 'i_2o keV — 
2.1 x 10 -9 erg cm -2 s _1 ) the source softens and the X-ray 
spectrum is better described by a —6 keV bremsstrahlung 
(Barret et al. 1998, in preparation). 

Terzan 2 is a relatively dim globular cluster, with 
an estimated distance of —10 kpc and visual extinction 
Ay = 4.03, according to Djorgovski (1993). The scatter on 
the values reported in the literature for these quantities is 
however rather high (~ 30% and ±5 5 %, respectively; see 
Ortolani et al. 1997 and references therein). The power- 
law intensity profile of its surface brightness and the rela- 
tively low core star density argue for a complex dynamical 
history (Djorgovski et al. 1986). 

Terzan 2 is one of the two globular clusters from which 
X-ray emission was detected at energies >40 keV by Sigma 
(Barret et al. 1991). The Sigma source (GRS 1724-308) 
was observed several times between 1991 and 1994 and 
appears to be variable (Goldwurm et al. 1994). The 1991 
detections (at a total 6.8er confidence level) are consistent 
with a simple power-law spectrum with T = 
38-200 keV flux of -5.8 x 10- 3 photons cm 
sponding to a poorly determined 1-20 keV flux of between 
1.5 and 13 x 10~ 10 erg cnr 2 s _1 . ROSAT HRI observa- 
tions by Mereghetti et al. (1995) did not reveal any other 
X-ray sources within 3'.5 from the cluster core, support- 
ing the identification of GRS 1724-308 with X 1724-308. 
The existence of a population a low- luminosity (i.e.: £ 
a few 10 33 erg s _1 ) sources in Terzan 2, similar to those 
discovered in the cores of 47 Tuc (Hasinger et al. 1994) 
or u> Cen (Johnston et al. 1994) cannot be excluded. The 
X-ray flux of these sources would be ~3 orders of mag- 
nitude lower than that of X 1724—308 and their existence 
therefore would not affect any of the conclusions in this 
paper. 

Recent RXTE observations of X 1724-308 (Olive et al. 
1998) suggest that its timing properties are typical of an 
"atoll" source. The Power Spectrum Density (PSD) can 
be well modeled with the sum of two "shot noise" com- 
ponents, with shot decay timescales of ~16 and ~680 ms. 
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A Quasi-Periodic Oscillation (QPO) feature at ^0.8 Hz 
was also detected in the PSD. The claim of another QPO 
feature at v ~ 0.092 Hz observed by the Medium Energy 
(ME) instrument on board EXOSAT (Belli et al. 1986) re- 
mains unconfirmed. The X-ray spectral hardness and the 
continuum properties of the PSD resemble those of some 
black hole candidates. 

As part of the BeppoSAX Core Program, a number 
of X-ray luminous globular clusters are being observed. A 
discussion of the properties of the sample will be deferred 
to a forthcoming paper. Here, we report the results of the 
BeppoSAX observation of X 1724-308 . 

2. Observation and data Reduction 

Results from all the co-aligned instruments on board Bep- 
poSAX (Boella et al. 1997a) are presented: the Low- 
Energy Concentrator Spectrometer (LECS; 0.1-10 keV; 
Parmar et al. 1997), the Medium-Energy Concentrator 
Spectrometer (MECS; 1.3-10 keV; Boella et al. 1997b), 
the High Pressure Gas Scintillation proportional Counter 
(HPGSPC, 4-120 keV, Manzo et al. 1997) and the 
Phoswich Detection System (PDS; 15-300 keV; Frontera 
et al. 1997). The MECS consisted at the time of the 
X 1724—308 observation of three identical grazing inci- 
dence telescopes with imaging gas scintillation propor- 
tional counters in their focal planes (now only two are op- 
crating) . The LECS uses an identical concentrator system 
as the MECS, but utilizes an ultra-thin (1.25 /im) entrance 
window and a driftlcss configuration to extend the low- 
energy response to 0.1 keV. The fields of view (FOV) of the 
LECS and MECS arc circular with diameters of 37' and 
56', respectively. The non-imaging PDS consists of four in- 
dependent units arranged in pairs each having a separate 
collimator. Each collimator can be alternatively rocked on- 
and off-source to monitor the background counting rate. 
The hexagonal PDS FOV is 78' full-width at half max- 
imum. The HPGSPC was similarly operating in rocking 
mode, but only the "off-" position has been employed to 
monitor the background. Several pieces of evidence indi- 
cate that there are no significant contaminating sources 
present in the PDS FOV, both in the "on" and in the 
"off" positions: the "off" count rates are in good agree- 
ment with the long-term increasing trend of instrumental 
background studied by Guainazzi & Matteuzzi (1997); the 
"on" and "off" light curves show the same variability pat- 
tern (~20% of count rates in the whole energy bandpass), 
which is due to the variation of the particle background 
around the orbit, thus excluding the presence of a strongly 
variable or transient unknown source and finally the PDS 
spectrum aligns very well with those from the other in- 
struments (see Sect. 3). Similar arguments apply to the 
HPGSPC data. 

The region of sky containing X 1724—308 was observed 
by BeppoSAX between 1996 August 17 04:29:05 UT and 
August 18 05:05:11 UTC in an observation lasting 65 ks 



and spanning 12 orbital cycles. Good data were selected 
from intervals when the elevation angle above the Earth's 
limb was >4° and when the instrument configurations 
were nominal, using the SAXDAS 1.3.0 data analysis 
package (Lammers 1997; custom software has been used to 
reduce the HPGSPC data). Additionally, PDS data within 
5 minutes after each South Atlantic Anomaly passage were 
excluded to avoid intervals were the gain was recovering to 
its nominal value following switch-on. This gives exposure 
times of 7.3 ks, 37.1 ks, 10.5 ks and 16.7 ks for the LECS, 
MECS, HPGSPC, and PDS, respectively. The standard 
PDS (HPGSPC) collimator rocking angle of 210' (180'), 
and standard dwell time of 96 s for each on- and off-source 
position were used. The LECS was operated only during 
satellite night-time and was switched-off during the second 
half of the observation for technical reasons. 



3. Data analysis 

LECS and MECS data were extracted centered on the po- 
sition of X 1724—308 using radii of 8' and 4' respectively, 
corresponding to ~ 95% of the instrumental Point Spread 
Functions. Background subtraction in the imaging instru- 
ments was performed using standard files, but is not crit- 
ical for such a bright source. Background subtraction in 
the HPGSPC and PDS was performed using data from the 
offset detectors. The background-subtracted count rates of 
X 1724-308 are 3.2 s"\ 10.4 s" 1 , 11.2 s" 1 and 6.1 s" 1 in 
the LECS, MECS, HPGSPC, and PDS, respectively. Spec- 
tra have been rebinned so as to have at least 30 counts per 
energy channel and to sample the instrumental energy res- 
olution with no more than 3 spectral channels, in order to 
ensure the applicability of x 2 fitting techniques. A sys- 
tematic error of 1% has been added to each channel of 
the rebinned LECS and MECS spectra, to account for the 
residual systematic uncertainties in the detector calibra- 
tion. Publicly available responses at 1997 September have 
been used, and spectral fits performed in the following 
energy ranges: 0.1-4 keV (LECS), 2.2-10.5 keV (MECS), 
4-40 keV (HPGSPC), 15-200 keV (PDS). In the spectral 
fits, a normalization factor has been included to account 
for the well-known mismatch in the BeppoSAX instrument 
absolute flux calibration (Cusumano et al. 1998). The fac- 
tors have been left completely free, except for the PDS 
to MECS one, which has been constrained to be in the 
0.82-1.02 range (Cusumano et al. 1998). The values of the 
LECS to MECS and HPGSPC to MECS relative normal- 
izations derived from the fit lie in the intervals 0.80-0.82 
and 1.05-1.06, respectively. They are in good agreement 
with the values typically observed (0.7-1.0 and 0.8-1.05, 
respectively; cf. Grandi et al. 1997; Cusumano et al. 1998). 
Uncertainties are given at 90% confidence level for one in- 
teresting parameter, unless otherwise specified. 
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Fig. 1. Light curves in the 0.1-1.5 keV (upper panel), 1.5- 
10.5 keV (central panel) and 13-200 keV (lower panel) en- 
ergy ranges. The binning time is 128 s. The LECS and 
MECS light curves are not background-subtracted. Only 
time bins with an exposure fraction > 25% are shown. 
Data from the two orbits after a Type I burst (correspond- 
ing to 12800 to 24200 s) have been removed 

3.1. Timing analysis of the persistent emission 

The source underwent a type I burst starting from 
~12825 s from the beginning of the observation, with an 
e-folding time of 25 ±3 s. In this paper we define the "per- 
sistent" emission excluding the data from the occurrence 
of the burst up to the end of the next orbit. The cut corre- 
sponds to ~10 4 s of elapsed time. Fig. shows light curves 
of the so defined persistent emission in the 1.5-10.5 keV 
(MECS), 0.1-1.5 keV (LECS) and 13-200 keV (PDS) en- 
ergy ranges. These energy ranges have been chosen in or- 
der to sample different spectral components (see Sect. 3.3). 
The x 2 °f a constant fit for these light curves is 245 for 
220 degrees of freedom (dof), 96.5/62 dof and 252/210 dof, 
respectively. At energies £ 1.5 keV there is no evidence of 
variability, while the relatively high x 2 obtained from the 
light curves at energy <1.5 keV is mainly due to an in- 
tensity increase of about 15% between the pre-burst and 
the post-burst phases. However, we will in the following 
analyze the time-averaged spectra of all instruments. In 
Sect 3.4 we show that the induced systematic uncertain- 
ties are smaller that the statistical errors on the spectral 
parameters, even in the simplest parameterization. 

3.2. The burst 

A profile of the burst in the full MECS bandpass is shown 
in Fig. |2|. The ratemeter CSELO (which registers the total 
count rate above a low energy threshold) is used to pre- 
serve the continuity of the burst profile in the plot, since 
some of the event scientific packets at the peak of the burst 
were missing due to telemetry losses. The burst shows a 
typical exponential decay with e-folding time 25 ± 3 s. We 



Fig. 2. Zoom of the MECS CSELO ratemeter light curve 
around the X1724-308 burst. Dashed vertical lines indi- 
cates the contiguous time intervals used for the time- 
resolved burst spectroscopy. For sake of clarity only a frac- 
tion of the "G" phase is shown, which extends for ~ 400 s 
from t ~ 13000 s 

have performed a spectral analysis of the burst in 7 con- 
tiguous time intervals (labeled as phases "A" to "G" in 
the following) . The borders of the intervals are shown as 
dashed vertical lines in Fig. |. The LECS and HPGSPC 
were switched off during most of the burst and PDS de- 
tection is generally too faint to allow time-resolved spec- 
troscopy on these small time scales. Therefore, spectra of 
the MECS data only have been extracted, except for the 
phases "F" and "G" , when LECS data were available. The 
average persistent spectrum has been used as background. 
Since the persistent emission in the 2-10 keV band is 
constant on timescales ~10 2 s within ±5%, time-resolved 
burst spectra were extracted under the condition that that 
2-10 keV count rate was £ 20% of the persistent one. At 
lower count rates, the uncertainties on the true instanta- 
neous level of the persistent emission make the determina- 
tion of the spectral parameters meaningless. The spectra 
have been fitted with a simple blackbody emission, photo- 
electrically absorbed by a column density, Nh, of neutral 
matter. All the parameters were allowed to vary, except 
the absorbing column density in the phase "C" spectral 
fit, which turned out to be unconstrained and was there- 
fore fixed to the average of the best-fit values in phases 
"E" and "F" and "G" . Best-fit parameters and results are 
reported in Table [l| The blackbody temperature decreases 
along the burst profile, as typically observed in this kind of 
burst. The apparent blackbody radius (obtained assuming 
spherical geometry) is consistent with being the same in 
all phases within the statistical uncertainties, except for 
the first ("A") and the latest ("F") phases, for which the 
best-fit is ~ 5 — 6 km. While the value obtained in the for- 
mer phase is likely to be strongly affected by the mixing of 
blackbodies with a rapidly changing distribution of tem- 
peratures during the sharp rise of the burst (and therefore 
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Table 1. Best-fit parameters and \ 2 values when a simple photoelectric absorbed blackbody model is applied to the 
time-resolved spectra of the X1724-308 burst. T exp is the exposure time, CR the (persistent emission subtracted) count 
rate in the 2-10 keV MECS bandpass, Rbb is the radius of a spherical blackbody emitting region and Lx the total 
X-ray luminosity. Errors on CR correspond to 1-a. 



Phase 


T CX p 


CR 


N H 


kT 


Rbb 


Lx 


x7dof 




(8) 




10 22 cm -2 


(keV) 


(km) 


(10 38 erg s _1 ) 




A 


2.8 


149 ±8 


< 6.3 


3.1±^ 


5± 2 , 


3-5±o;7 


4.6/9 


B 


1.17 


420 ± 20 


< 4.0 


2.2±8:i 


13±| 


5.4± f 1 , : ? 


6.6/10 


C 


0.79 


480 ± 20 


1.2* 


i.93±8:?g 


18 ±4 


6.0 ±0.2 


5.5/8 


D 


3.9 


250 ±8 


< 1.7 


2.i±8:i 


11 ±3 


3.i±8:l 


25/20 


E 


31.7 


78.0 ±1.7 


1.4 ±0.5 


1.36 ±0.06 


12.5±l;i 


0.69 ±0.04 


41.2/40 


F 


104.7 


18.8 ±0.5 


i.3±8:? 


1.08 ±0.08 


9.2±ii 


o.i5i±8:8ii 


60.3/40 


G 


401.9 


3.22 ±0.19 


o.8±8:l 


0.87±S:i ( 1 , 


5.8±U 


0.025 ±0.004 


65.9/65 



' fixed 



probably unreliable), the latter could suggest a true de- 
crease of the emitting region typical size with time; this 
suggestion is supported by the fact that the nominal best- 
fit value in phase "F" is indeed lower than between phases 
"B" to "E" . However, better statistical quality of the data 
is required to confirm such a conclusion. The weighted av- 
erage of the Rbb is < Rbb >= 9.9 ± 0.7 km (11.6 ± 0.8 km 
if only phases between "B" and "F" are considered). The 
peak flux is ~ 5.0 x 10~ 8 erg cm~ 2 s , corresponding to 
a luminosity ~ 6 x 10 38 erg s . The simple parameteri- 
zation adopted yields an adequate fit for all spectra. 




3. 3. Spectral analysis of the persistent emission 
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We performed the spectral analysis of the persis- 
tent emission in several steps. First LECS/MECS and 
HPGSPC/PDS spectra were fitted separately with purely 
phcnomenological models. The spectra of all detectors 
were then fitted simultaneously, assuming different physi- 
cal scenarios. 

A single component model with absorption provides 
a rather poor description of the LECS/MECS spectrum. 
Power-law, bremsstrahlung and blackbody spectral mod- 
els were tried. All of these give large excesses above the 
extrapolations of the high-energy data at ~1 keV (see 
Fig. |^) . A two-component model is needed to account for 
the observed spectrum (see Table ^). Formally, the best 
combination is a relatively flat (r ~ 1.6) power-law ± a 
blackbody with kT ~ 1 keV. Other combinations yield a 
null hypothesis likelihood ~ 1%. 

X 1724-308 is detected in the PDS up to ~ 150 keV. 
One-component models give similarly poor descriptions 
of the HPGSPC and PDS spectra (see Table |). There is 
strong evidence of a sharp roll over of the PDS spectrum 
above ~ 40 keV (see Fig. ^) . The inclusion of a cutoff in the 
power-law model significantly improves the quality of the 
fit (Ax 2 = 129, xl — 1-06), yielding best-fit parameters of 
r = 1.86±0.04 and E cutoff = 90±}§ keV. There is evidence 
in the HPGSPC data that the spectral index undergoes a 



Fig. 3. LECS and MECS spectra (upper panel) and 
data/model ratio (lower panel) when a photoelectric ab- 
sorbed power-law is applied 




50 

Energy (keV) 



Fig. 4. Ratio of the PDS spectra of Terzan 2 and Crab 
(upper panel) and difference residuals (lower panel) when 
a power-law model is applied to fit the ratio in the 13- 
30 keV band and extrapolated to higher energies (solid 
line). Each data point has a S/N ratio > 20 



Evidence of Comptonization from the broadband spectrum and variability of the X-ray source in Terzan 2 



5 



Table 2. Best-fit parameters and \ 2 values when the models in the first column are applied simultaneously to the LECS 
and MECS spectra. In column 1 or in the subscripts: wa = photoelectric absorption; po — power-law; bb — blackbody; 
brems = thermal bremsstrahlung. 



Model 


N H 


r Or kTbrems 


Tbrems Or Tbb Or Ebreak 


X7 dof 




(10 22 cm" 2 ) 


(keV) 


(keV) 




wa*po 


1.86 ±0.05 


1.834±olooi3 




303.4/76 


wa*bknpo 


i.42±8:8? 


1 e:c:_i_0.04 
1.00=t(j.o6 

i-9i±8:81 


4.3±8:1 


108.4/74 


wa* (po+bb) 


1.31 ±0.08 


1.61 ±0.05 


i.08±8:8S 


89.5/74 


wa* (po+brems) 


2.39 ±0.09 


1.924 ±0.018 


0.18 ±0.02 


101.6/74 


wa*brems 


1.38 ± 0.04 


n.6±8:l 




123.5/76 


wa* (brems+bb) 


l-46± .05 


n.2±8:i 


0.10 ±0.03 


104.0/74 



Table 3. Best- fit parameters and \ 2 values when the models in the first column are applied simultaneously to the 
HPGSPC and PDS spectra, cutoffpl = power-law± cutoff. 



Model 


r 


kTbr 

(keV) 


Ecutoff 

(keV) 


X7 dof 


po 


2.104 ±0.014 






266.6/131 


brems 




28.2±g;l 




667.6/131 


cutoffpl 


1.86 ±0.04 




99±ii 


137.6/130 


po+brems 


2.14 ±0.08 


30±| 




171.7/129 



trum: an intermediate-hard (i.e.: 1.5-50 keV) power-law 
"emission core", with spectral index ~ 1.6 — 1.9 (proba- 
bly gently steepening with energy); a sharp cutoff above 
~ 50 keV and a soft excess at energies below ~ 1 keV. 

The models proposed so far to explain the X-ray emis- 
sion from low-mass X-ray binaries (LMXRB) can be cate- 
gorized in two main classes. The first one (originally pro- 
posed by Mitsuda et al. 1984) postulates the superposi- 
tion of a multicolor blackbody emission from the inner 
region of an accretion disk (typical effective temperature 
Tdisk ~1 keV) and of a simple blackbody originating on 
the surface of the neutron star (Tbb ~2 keV). Another 
class invokes unsaturated Comptonization of low-energy 
photons (White ct al. 1985, 1986), possibly coupled with 
a blackbody (White et al. 1988). It has been recently sug- 
gested that these models characterize different luminosity 
states in LMXRB (see e.g. Tanaka 1995), the former pre- 
vailing for Lx ~ 10 37 erg s _1 . Actually, the double black- 
body scenario is ruled out by the present data, at least 
in its simplest form. The fit with such a model yields a 
very poor \ 2 (xi = 21); the proposed model does not 
produce enough flux at high energy to account for the 
relative hardness of the observed spectrum. 

We have tested the viability of the Comptonization 
scenario with the xspec model compTT (Titarchuk 1994; 
Hua & Titarchuk 1995; Titarchuk & Lyubarskij 1995), 
which calculates self-consistently the spectrum produced 
by the Comptonization of soft photons in a hot plasma. 
This model contains as free parameters the temperature of 
the Comptonized electrons T c , the plasma optical depth r 
and the input temperature of the soft photon (Wien) dis- 
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Fig. 5. In the upper panel the ratio of the HPGSPC spec- 
tra of XI 724-308 and the Crab Nebula multiplied by 
E~ 21 (typical value of Crab Nebula photon index, see 
e.g. Cusumano et al., 1998). The solid line is the best-fit 
with the energy dependency law: 0.19 x J5-l-72-o.ulo g (s) j 
The residuals in ratio between the normalized Crab ra- 
tio spectrum and the best-fit law above are shown in the 
lower panel 



smooth increase with energy (see Fig. ||) from T ~ 1.80 at 
5 keV to T ~ 1.90 at 30 keV. 



3.3.1. Comparing physical scenarios 

Three spectral region can be phenomenologically recog- 
nized in the X 1724— 308 broad-band BeppoSAX spec- 
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Table 4. Best-fit results when Comptonization models are applied simultaneously to the LECS, MECS and PDS 
spectra, diskbb = multi-temperature accretion disk compTT = Comptonization. For the compTT model the letter in 
bracket indicates the disk ("d") or spherical ("s") geometry for the Comptonizing plasma, y is the Comptonization 
parameter (= 4[kT c /m c c 2 ]r 2 ]) and Rw is the typical size of the input photons emitting region (following Titarchuk 
1994) 



# 


Model 


N H 


kT c 


T 


y 


kT w 


Rw 


kT sc 


R 


x7dof 






(10 22 cm" 2 ) 


(keV) 






(keV) 


(km) 


(keV) 


(km) 




1 


wa* (compTT(d)) 


2.06±8:ffi 


24±| 


i.42±&i? 


0.4 


0.15±H 2 


350 






518.6/203 


2 


wa* (compTT(d)+bb) 


o.78±8:J? 


29±3 8 


i.i7±8:S 


0.3 


i.02±8il 


7.9 


0.60±8;8| 


12±| 


199.1/201 


3 


wa* (compTT(d)+diskbb) 


i.n±8:8S 


28±s 3 


i-3±8;2 


0.4 


2.io±8i| 


1.8 


i.5o±8:?l 


2.5±g : r 


194.4/201 


4 


wa*(compTT(s)) 


0.48±8:8a 


23±t 


3.7±8:2 


2.5 


o.720±8:8S8 


15 






337.6/203 


5 


wa*(compTT(s)+bb) 


o.78±8:8? 


30±1 


2.9±8:S 


2.0 


i.oa±8:io 


8.3 


o.6o±8:81 


12 ±3 


198.5/201 


6 


wa* (compTT(s)+diskbb) 


1.13 ±0.04 


27±i 2 


3.3±8:f 


2.3 


i.98±8:?g 


2.1 


i-44±8:Si 


2.7±8:f 


194.9/201 



R = R' cos(9) 1 ^ 2 , where 9 is the inclination angle and R' is the inner radius of the disk 
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Fig. 7. Left: LECS, MECS, HPGSPC and PDS spectra of X 1724-308 (upper panel) and residuals in units of standard 
deviations (lower panel) when a the best-fit model #5 in Table 3 is applied. Each data points has a S/N ratio > 3. 
Right: the inferred photon spectrum with the best-fit model, deconvolved in the Comptonized (dashed line) and disk 
(dot- dashed line) blackbody components 



tribution Tyj. Both a disk and a spherical geometry for 
the Comptonizing plasma were tried. Again, a soft excess 
remains around ~ 1 keV if compTT alone is applied. We 
have tried to model this excess either with a single tem- 
perature or with a multi temperature blackbody, for each 
of the two plasma geometries. The soft excess cannot in- 
stead be modeled with a Gaussian line, even if its intrinsic 
width is allowed to be broad. The results are summarized 
in Table ||. The above scenarios cannot be discriminated 
on purely statistical grounds, since the values of xt are 
all very close to one. The Comptonized plasma is char- 
acterized by an electron temperature ~25-35 keV (90% 
confidence level for two interesting parameters, see Fig. ^]) 
and a geometry- independent parameters (3 ~ 0.2, which 
expresses the probability of photon scattering through the 
exponential factor e _/3 (Titarchuk 1994). The other fit pa- 
rameters are slightly dependent on the assumed model 
used to describe the soft X-ray excess, typical values be- 



ing N H ?t 0.8 (1.1) x 10 22 cm- 2 , kT w ~ 
kT S0 ~ 0.6 (1.5) keV and R ~ 12 (3) km, 



1 (2) keV 
when it 



is 




25 



30 
kT (keV) 



35 



Fig. 6. Contour plot r vs. kT e for the model #5 in Table |[ 
Iso-x 2 curves are at 90% and 99% level of confidence for 
two interesting parameters 



described with a single- (multi-) temperature blackbody. 
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The best-fit model #5 in Table |] superposed to the ob- 
served spectra is shown in Fig. |7|. It corresponds to the fol- 
lowing fluxes: 6.4 x ICT 10 erg cm" 2 s" 1 (2-10 keV), 10.3 x 
10-10 erg cm -2 s -i ^_ 20 keV ) 7 4.7 x lO" 10 erg cm" 2 s" 1 

(38-200 keV), 18.3 x 10~ 10 erg cm- 2 s" 1 (0.1-200 keV), 
the last one corresponding to an unabsorbed X-ray lumi- 
nosity Lx ~ 2.3 x 10 37 erg s^ 1 for a distance of 10 kpc 
(Djorgovski 1993). Given such a continuum, the 90% up- 
per limit on the equivalent width of a narrow line at 
6.7 keV is 17 eV. 

However, it must be noted that the broadband spec- 
trum can be formally accounted by a suitable combina- 
tion of three thermal models, as shown in Table |5|. A 
bremsstrahlung with temperature Tb r — 45 — 60 keV is 
required to fit the PDS hard tail, a 0.8-1 keV blackbody 
can account for the soft excess, while a 1.9 keV blackbody 
or a ~ 9 keV bremsstrahlung are needed to fill the gap 
between the two "extreme" spectral components. In both 
cases, the x 2 is as good as for the Comptonization models. 

3.3.2. The soft X-ray variability 

After the burst, the 0.1-1.5 keV flux increased by ~ 15%. 
LECS and MECS spectra corresponding to these inter- 
vals were extracted and analyzed with the same model 
(model #5 in Table | with kT c = 28 keV and r = 3.1), 
in order to investigate if there is any significant spectral 
variability associated with this flux change. The best-fit 
parameters are summarized in Table ^[ The x 2 f°r t ne 
combined fit is acceptable (x 2 = 1.13). All the best-fit 
parameters are mutually consistent between the two in- 
tervals. We conclude that the systematic uncertainties in- 
troduced by mixing models with different normalizations 
are well below the statistical uncertainties, given the qual- 
ity of the current data. 

4. Discussion 

We report on the first X-ray broadband spectrum of the 
burster X1724-308 in the whole 0.1-100 keV domain. The 
source is detected with a S/N > 3 up to ~150 keV, con- 
firming the previous Sigma/GRANAT results (Barret et 
al. 1991). The superior combination of broadband cov- 
erage, sensitivity and energy resolution has however al- 
lowed an unprecedentedly detailed spectral deconvolution, 
albeit formally not unique. Qualitatively the spectrum 
can be described as the composition of a "core", power- 
law emission, which undergoes a significant steepening at 
E~ 50 keV, and a soft excess above the extrapolation of 
the higher energy X-ray spectrum. The 20-200 keV flux is 
~ 3.4 x 10~ 10 erg cm~ 2 s _1 . The corresponding luminos- 
ity at 10 kpc is ~ 4.1 x 10 36 erg s _1 , which is typical of 
X-ray bursters (van Paradijs & van der Klis 1994; Barret 
et al. 1996). 

Formally, the spectrum can be described by a suit- 
able combination of three thermal models (models #7 



and #8 in Table The characterizing feature of these 
models is the need for a very high temperature thermal 
bremsstrahlung (Tb r ~ 50 — 60 keV) to account for the 
highest energy PDS data points and cut-off roll over. The 
presence of such a component had already been invoked 
to explain the hard tail observed in several X-ray bursters 
(Churazov et al. 1995). The luminosity inferred by the 
BeppoSAX observation is Lb r ~ 1.4 x 10 37 erg s _1 , which 
implies an emission measure < n V 2 > ~ 3 x 10 59 cm 3 ; 
the optically thin requirement leads to a size of a single 
uniform emitting cloud >7 x 10 10 cm. Such a large region 
is unlikely to be directly powered by the neutron star. 

The high-energy tail rules out the "double-blackbody" 
model (multi-temperature disk + neutron star surface), 
which had originally proposed by Mitsuda et al. (1984) 
to explain the X-ray emission of LMXRB. If a transi- 
tion from a thermal- to a Comptonization-dominated hard 
X-ray spectrum is present in X 1724— 308 (as generally 
believed for LMXRB), it occurs at luminosity £ a few 
10 37 erg s _1 . 

Models of the BeppoSAX spectrum of X 1724-308 
which involve a Comptonized component appear there- 
fore more plausible. A plasma with kT e ^27-30 keV and 
(3 ~ 0.2 is required. A disk and a spherical geometry of 
the Comptonizing region fit the data comparably well. 
We note however, that somewhat higher values of the 
Comptonization parameter (y ~2.0-2.3) are obtained for 
a spherical geometry. Similar Comptonized spectral com- 
ponents have been observed in a variety of neutron star 
LMXRBs, although in most cases the spectra are substan- 
tially softer and the inferred electron temperatures in the 
2-4 keV range (White et al. 1988). In a few Atoll sources, 
however, X-ray spectra extending up to ^100 keV or more 
have been revealed through Sigma/GRANAT (Mandrou 
et al. 1994) and BeppoSAX (Church et al. 1998b) obser- 
vations. These spectra have a clear resemblance to the 
spectrum of X 1724—308 discussed here, and require com- 
parably high electron temperatures. It has also to be con- 
sidered that hard tails extending to hundreds of keV are no 
longer believed to be exclusive to black hole candidates, 
at least when the LMXRB luminosity does not exceeds 
~10 37 erg s _1 (Barret et al. 1994; van Paradijs & van der 
Klis 1994). 

In the Beppo-SAX observation of X 1724-308 an ad- 
ditional soft spectral component is needed to fit the low 
energy spectrum. However, the quality of the data does 
not allow different models for this component to be dis- 
criminated. We find that a multi-temperature blackbody 
disk model, or a simple blackbody spectrum provide ac- 
ceptable fits. 

In the case of a multi-temperature blackbody disk 
model, both the derived temperature (~2 keV) and ra- 
dius of the innermost disk region (~3/(cos0) 1 / 2 km) are 
close to those obtained for the "seed" Wien spectrum for 
the Comptonized component (see e.g. the model #6 in 
Table |J) . An inclination angle of ~ 70° would be required 
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Table 5. Best-fit results for the "thermal scenarios" of the Terzan 2 broadband spectrum. The subscripts of the kT 
refer to the thermal components in the model in column 1 read from the left to the right. Rbb is the size of a spherical 
optically thick blackbody emitting region 



# 


Model 


N H 
(10 22 cm" 2 ) 


(keV) 


Rbb 

(km) 


kT 2 
(keV) 


kT 3 
(keV) 


X Vdof 


7 

8 


wa* (bb+bb+brems) 
wa* (bb+brems+brems) 


o.95±g:S 
i.25±8:8S 


o.83±8:£ 
o.99±8:8S 


6.5±^ 

i.i7±8:i8 
2.s±8:t 


i.98±8:i| 
9.3±i:§ 


48 ±3 
62±? 


212.3/201 
202.2/201 



Table 6. Best-fit parameters for the joint fit of the LECS and MECS spectra extracted before and after the burst 



Temporal phase 


N H 


kT w 


kTbb 


N bb 




(10 22 cm" 2 ) 


(keV) 


(keV) 


(10~ 3 L 39orgs -iD M 2 kpc ) 


Before the burst 


087±;j:^ 


1.1 ± 0.2 


0.64±!J;i!l 


2.6±H 


After the burst 


0.84±8;J| 


0.99±»;}? 


0.57±°;° 7 


2.3 ±0.5 



for the inner radius of the disk to accommodate a neutron 
star. On the contrary the absence of a sizeable orbital 
modulation of the X-ray flux (absorption dips in particu- 
lar) suggests a substantially lower inclination (P89; Olive 
et al. 1998). 

If the soft thermal emission of Terzan 2 is interpreted in 
terms of a single temperature blackbody, the inferred lin- 
ear size of the emitting region (~ 12 km) is in good agree- 
ment with a typical neutron star and with the size inferred 
from the analysis of the type I burst occurred during the 
BeppoSAX observation. It is worth emphasizing the (sin- 
gle) blackbody temperature of Terzan 2 (kT~ 0.6 keV) is 
a factor of 2-3 lower than that derived from most other 
neutron LMXRBs (cf. White et al. 1988). In the single 
temperature blackbody scenario, the parameters also are 
close to those derived for the input spectrum to the Comp- 
tonizcd component, such that the neutron stellar surface 
might represent the source of "seed" photons. Alterna- 
tively, the analogy with the hard spectra of black hole 
candidates (where, of course, there is no star surface or 
boundary layer) suggest that even in disks around accret- 
ing neutron stars there is a hot Comptonizing plasma that 
coexists with a dense cold phase responsible for the gener- 
ation of the "seed" photons (cf. Eardley & Lightman 1975; 
Galeev et al. 1979; Molcndi & Maraschi 1990). 

The two spectral decompositions discussed above dif- 
fer a great deal in the physical interpretation of the two 
spectral components. In the former case the Comptonizcd 
hard component originates from the boundary layer close 
to the neutron star surface, whereas the soft component 
is produced by the sum of blackbody spectra emitted by 
the accretion disk. In the latter case the (single) black- 
body spectrum likely arises from the boundary layer close 
to the neutron star surface, while the accretion disk gives 
rise to the Comptonized spectrum. Simple accretion disk 
theory predicts that, unless the neutron star is rotating 
very close to its break up velocity, the ratio of the accre- 



tion luminosity released in the disk to that released in the 
boundary layer should be equal to 1. Sunyaev & Shakura 
(1986) show that this ratio can be as low as ~ 0.45 if 
the accretion disk ends at the marginally stable orbit and 
this, in turn, is larger then the neutron star radius. Re- 
cent observations of fast quasi periodic oscillations, QPOs, 
in LMXRBs testify that the neutron stars contained in 
them rotate much slower than break up (periods in the 
—3 ms range), and that the inner radius of the accre- 
tion disk (as inferred from the frequency of the higher 
kHz QPO power spectrum peak) is often larger than the 
marginally stable orbit. This would further decrease the 
predicted ratio of disk and boundary layer luminosity. The 
0.1-100 keV luminosities of the disk (multi-) blackbody 
component and the (boundary layer) Comptonized com- 
ponent inferred from the BeppoSAX spectrum of Terzan 
2 are -8.5 x 10 36 erg s" 1 and - 1.5 x 10 37 erg s" 1 . The 
derived ratio of —0.55 is in agreement with the predictions 
above. However, observations of dips in several X-ray can 
be convincingly explained if the Comptonized component 
is extended (Church & Balucihska-Church 1995; Church et 
al. 1998a), almost at odds with the former interpretation. 

In the latter interpretation the ratio of the luminosity 
of the Comptonized (disk) spectrum (1.9 x 10 37 erg s _1 ) 
and the luminosity of the (boundary layer) single temper- 
ature blackbody (2.4 x 10 36 erg s _1 ) is — 8, i.e. much 
higher then expected. An anomalously high ratio of the 
two components has also been found in a number of other 
neutron star LMXRBs (see White et al. 1988). Obscu- 
ration by an inner accretion corona or a thickened disk 
could decrease the relative strength of the boundary layer 
contribution (Lamb 1986; van der Klis et al. 1987). In 
the standard Shakura-Sunyaev disk, the inner radiation 
pressure-dominated disk can be as thick as ~ few tens of 
kilometers, and the effect of thermal and/or secular insta- 
bilities can further enlarge it (Shakura & Sunyaev 1973; 
Pringle 1981). Radiation from the boundary layer might 
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Fig. 8. Beppo-SAX and EXOSAT count spectra. EX- 
OSAT data point have been multiplied by 2.5 to align 
the spectra below the photoelectric absorption cut-off 

be partly scattered out of the polar channel or attenuated 
through a not completely optically-thick disk. A "bloated- 
disk" geometry has already been proposed to explain the 
low frequency QPOs observed in several LMXRBs (e.g.: 
4U1820-30, Stella et al. 1987), and would be in line with 
the likely low inclination of the system, as inferred from 
the lack of significant intensity eclipses. Alternatively, ex- 
change of energy between the magnetosphere and the ac- 
cretion disk via material torque from the accretion matter 
could reduce the luminosity available at the surface of a 
neutron star, spinning at a rate close to the inner disk 
Keplerian frequency (Priedhorsky 1986). 

The optical reddening is Ay = 4.03 (Djorgovski 1993). 
If we convert the visual absorption in an X-ray absorb- 
ing column density with Nh = 1.8 x 10 21 Ay cm -2 (Pre- 
dehl & Schmitt 1995), the corresponding column density 
is roughly consistent with the Galactic one along the line 
of sight to Terzan 2 (N HGaJ = 6.5 x 10 21 cm" 2 , Dickey & 
Lockman 1990). However, the BeppoSAX observation re- 
quires a 20-75% excess absorption (depending on the spec- 
tral descriptions), in agreement with previous estimates, 
although rather poorly constrained and model dependent 
(P89; Mereghetti et al. 1995; Verbunt et al. 1995). 

EXOSAT observed X 1724-308 in a factor -5 fainter 
and in a much softer state than BeppoSAX. This is shown 
in Fig. ||, where the photon spectra observed by Bep- 
poSAX and EXOSAT arc compared. The hardness dif- 
ference can be parameterized as a Ar ~ 0.7 between the 
observed power law photon indices. If the electron tem- 
perature remains the same, (3 ~ 0.3 and therefore a ^50% 
smaller r (Hua & Titarchuk 1995) at the EXOSAT epoch 
is implied (we note however that in this regime the analyt- 
ical approximations on which the Comptonization Monte- 



Carlo simulation by Titarchuk (1994) are no longer valid 
and therefore the condition that T c and [3 determines uni- 
vocally the spectral index is not fulfilled, Hua & Titarchuk 
1995). Alternatively, a steeper spectrum could imply a 
lower kT c and therefore a spectral cut-off at lower en- 
ergy. To check this possibility, we extracted from the pub- 
lic archive and re-analyzed the EXOSAT ME spectrum of 
P89. A simple absorbed power-law model gives a reason- 
able fit (xt — 77/61 dof), with best-fit parameters as in 
P89. However, the inclusion of a cut-off results in a reduc- 
tion in x 2 of 8 which is significant at 98.9% confidence. The 
best-fit parameters are then N# = (1.5 ± 0.2) x 10 22 cm~ 2 
(very similar to the BeppoSAX value in model #6 here), 
r = 2.0 ± 0.2 and E cut _ off = 14± 21 keV. We caution that 
this cut-off energy is very close to the upper energy range 
of the ME detector. If the EXOSAT result is correct, then 
an explanation of the difference in the two spectra may be 
a factor ~ 4 lower electron temperature during the EX- 
OSAT epoch. In such a scenario the energy of the electrons 
is physically linked with the high-energy intensity. 
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